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The phase diagram of correlated, disordered electron systems is calculated within dynamical 
mean-field theory for the Anderson-Falicov-Kimball model with nearest-neighbors and next-nearest- 
neighbors hopping. The half-filled band is analyzed in terms of the chemical potential of the system 
using the geometric and arithmetic averages. We also introduce the on-site energies exhibiting 
a long-range correlated disorder, which generates a system with similar characteristics as the one 
created by a random independent variable distribution. A decrease in the correlated disorder reduces 
the extended phase. 



PACS numbers: 71.10.Fd, 71.27.+a, 71.30.+h 



I. INTRODUCTION 

Many recent works [l], [H, Q have considered a non- 
perturbative framework to investigate the Mott metal- 
insulator transition (MIT) Q in lattice electrons with lo- 
cal interaction and disorder using the dynamical mean- 
field theory (DMFT) @. For example, the Anderson 
transition jy, @] has been explored on the Bethe lattice 
Q considering the Hubbard and Falikov-Kimball models. 
In these studies, the Mott MIT is characterized by open- 
ing a gap in the density of states at the Fermi level. At 
the Anderson localization, the character of the spectrum 
at the Fermi level changes from a continuous to a dense 
discrete one. 

The study of disordered systems requires the use of 
probability distribution functions (PDFs). One is usu- 
ally interested in typical values of these quantities which 
are mathematically given by the most probable value of 
the PDF. The metal and the insulator phases could be 
detected by analyzing the local density of states (LDOS). 
In particular, the arithmetic mean of this random one- 
particle quantity is noncritical at the Anderson transition 
and hence cannot help to detect the localization transi- 
tion. By contrast, the geometric mean gives a better 
approximation of the averaged value of the LDOS [1, [t| , 
as it vanishes at a critical strength of the disorder and 
hence provides an explicit criterion for Anderson local- 
ization H, E3, Ellll- Recently, we adopted the Holder 
mean to analyze how the averaged LDOS depends on 
each Holder parameter that is used. We showed that the 
averaged LDOS can vanish in the band center at a criti- 
cal strength of the disorder for a wide variety of averages 

Most of these studies are restricted to the hopping be- 
tween nearest- neighbors (NNs). Recent works suggest 
that the inclusion of next-nearest-neighbors (NNNs) hop- 
ping in the model favors the long-range magnetic order- 
ing [14| . The investigation of these effect is particularly 
interesting. The breaking of particle-hole symmetry even 
at half filling, is a generic property of real materials and 



creates effects on the paramagnetic phase [15l|. b esides 
the frustration of the antiferromagnetic phase [16fl . The 
effect of nonrandom NNN hopping becomes evident al- 
ready in the noninteracting system through an asymmet- 
ric density of states (DOS), as alread y d erived for an ar- 
bitrary hopping on the Bethe lattice \u\ . 

On the other hand, the long-range correlated disor- 
der can be generated in a variety of stochastic processes 
in nature [18j |. A tight-binding one-dimensional model 
of electronic states with the on-site energies exhibiting 
long-range correlated disorder and nonrandom hopping 
amplitudes was studied in ref. [l9j |. The presence of an 
Anderson-like metal insulator transition was revealed for 
a finite range of energy values where the Lyapunov co- 
efficient vanishes. The correlation in the disorder favors 
the emergence of the extended phase. 

In this paper, we investigate two different aspects for 
the Anderson-Falicov-Kimball model. First, consider- 
ing nearest-neighbor (NNs) and next-nearest-neighbor 
(NNNs) hoppings, we analyzed how the presence of the 
NNN hoppings influence the phase diagram of the ground 
state of this model using geometric and arithmetic aver- 
ages. Besides that, we showed how the chemical potential 
varies as we look at the half-filled band of different sys- 
tems; the inclusion of the NNN hopping dislocates the 
half-filled band from the band center, independent of the 
disorder that is considered. The ground-state phase dia- 
gram is also presented for different values of NNN hop- 
pings. 

Secondly, we studied the main effects of the long-range 
correlated disorder (characterized by the exponent a) in 
the Anderson-Falicov-Kimball model. Unlike the one- 
dimensional Anderson model, this disorder does not have 
a great influence on the behavior of the system. We 
present the ground-state phase diagram for electrons in 
a half-filled band for different values of this disorder and 
the dependence of the disorder in terms of the exponent 
a for a system without Coulomb repulsion. 

In the pure half-filled Falicov-Kimball model, the 
Fermi energy for electrons is inside of the correlation 



(Mott) gap opened by increasing the interaction 0. We 
want to examine how the disorder influences this gap. 

The pure Falicov-Kimball model describes two species 
of particles, mobile and immobile, which interact with 
each other when both are on the same lattice site. This 
model is the simplest model to study metal-insulator 
transitions in mixed valence compounds of rare earth 
and transition metal oxides, ordering in mixed valence 
systems, order-disorder transitions in binary alloys, itin- 
erant magnetism [2fj, crystallization, electronic ferroelec- 
tricity in mixed- valence compounds, and phase diagrams 
of metal ammonia solutions [2l|. It also captures some 
aspects of the Mott-Hubbard MIT [Hj]. 

The Anderson- Falicov-Kimball model considers mobile 
particles that are disturbed by a local random potential, 
giving rise to a competition between interaction and dis- 
order. The Hamiltonian is written as 
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where (c,) and (fa) are, respectively, the cre- 
ation (annihilation) operators for the mobile and immo- 
bile fermions (electrons and ions, respectively) at a lattice 
site i, tij is the electron transfer integral connecting sites 
i and j, and U is the Coulomb repulsion that operates 
when one ion and one electron occupy the same site. We 
assume 

t\ for nearest-neighbor sites 
t2 for next-nearest-neighbor sites 
otherwise 

The average number of electrons and ions on site i are 
denoted, respectively, as n e = cfci and nj = /,•/». We 
consider that the occupation n/ on the zth site has prob- 
ability p (0 < p < 1). It was assumed, for simplicity, 
that just mobile particles are subjected to the structural 
disorder Q. The energy ej is a random, independent 
variable, describing the local disorder disturbing the mo- 
tion of electrons. The model is solved within the DMFT 
framework. 

In section [II] we present the DMFT approach applied 
to the Anderson- Falikov-Kimball model |9j . In section HTll 
we present the numerical results concerning the ground- 
state phase diagram for the next-nearest-neighbor hop- 
ping case. In section [IV] we discuss the above model for 
the long-range correlated disorder. Finally in section fVl 
we present our conclusions. 



II. DYNAMICAL MEAN-FIELD THEORY 

The formalism is based on the introduction of the hy- 
bridization function rj(uj), which is a dynamical mean 
field describing the coupling of a selected lattice site with 
the rest of the system [12J . The DMFT is calculated from 
the Hilbert transforms 
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FIG. 1: Averaged local density of states at U — 3.6 and q = 
for (a) tl/t* = 0, 0.4 and 1.0 for disorder strength A = 2.0. 



and 



(3) 



where Nq(c) is the non-interacting density of states and 
G(u>) the translationally invariant Green function. 
The e^-dependent LDOS is written as 



p(uj,e l ) = -- Im G(oj,€i 

IT 



(4) 



where G(u>, is the local e^-dependent Green func- 
tion Q. £i is considered an independent random vari- 
able characterized by a probability function P{ti) = 
$(A/2 — £i)/A, with <E> being the step function. The pa- 
rameter A is a measure for the disorder strength. From 
the ei-dependent LDOS, we introduce the (/-Holder aver- 
aged LDOS 
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where the subscript q defines the generalized mean. Spe- 
cial cases are, for example, the minimum {q — > -co), the 
geometric mean (q — ► 0), the arithmetic mean (q = 1) 
and the maximum (g — + oo). Eq. ([5]) inserted in Eq. ([3]) 
closes the self-consistent DMFT. 

A chemical potential /i is introduced only for the mo- 
bile subsystem, to fix the system in the half-filled band 
(n e = p = 1/2). For the Bcthc lattice in the limit of infi- 
nite connectivity K, one can use the scaling t\ — t*/^/K 
and ti = t* 2 /K. In the nearest-neighbor hopping case 
(<2 = 0) the density of states is the semi-elliptic function 
N„(e) = a/4 - (e/q) 2 /(27rfj) [H, where the bandwidth 
is 4t*. We can find r?(o>) = tfG(u>) 0. In this case, the 
ground-state properties in the half-filled band case are 
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FIG. 2: Phase diagram of the ground state for the Anderson- 
Falikov-Kimball model for U = 3.6 at = (solid line), 
iSA* = 0.4(dash-dot line) and ij/t* = 1 (dash-dot-dot line). 
Black lines present mobility edges determined within DMFT 
with geometric averaging (g = 0) and brown lines show band 
edges determined within DMFT with arithmetic averaging 

(9 = 1). 



solely determined by the quantum states in the band cen- 
ter (uj = 0), and we can determine the transition points 
in the phase diagram by linearizing the DMFT equations. 
However, for t 2 ^ the symmetry in the band center is 
absent, and we can not use a recursive relation within 
the linearized DMFT. 

For the next-nearest-neighbor hopping case (f 2 7^ 0) 
the analytical expression for the density of states on the 
Bethe lattice is given by 



N (e) = <d{tf + Atf + At* 2 e) 



A 2 



4t* 2 e 
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where 9(e) is the step function and A±(e) is solution of 
the equation 



t* 2 X 2 +t\\- (t* 2 +e) = 0. 



(7) 



III. NEXT-NEAREST-NEIGHBOR HOPPING 

We considered that the initial value of p(uj, e$) is a uni- 
form distribution with bandwidth W = 8i*, then we de- 
termined G(<jj) in order to obtain ?y(w) and finally the new 
values of p(u>,ei). This loop is performed until we find 
the stable configuration for p q (w). The relation between 
G(lo) and rj(oj) is obtained in a straightforward way from 
Eqs. ((2]) and ([6|) leading to the result 



r,(cj) = G(u) 



tf t?(l-G(m)t* 2 ) 
1 - G{u)t* 2 (1- 2G{io)t* 2 ) 2 
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If we consider fixed values of U and t 2 /t\ we can ob- 
serve the dependence of the averaged LDOS on the value 
of A that is used. Some of these results are presented in 
Fig. [TJ We consider the next-nearest-neighbor hopping 
model (t 2 /t* = 0.0, 0.4 and 1.0). We can see that the 
behavior of the LDOS is not the same in all cases. In 
contrast to the case t* 2 jt\ = 0, if t 2 /t* 7^ the p q {oj) is 
asymmetric. 

A signature of the Anderson localization is the van- 
ishing of p q {uj) as we increase A. The detection of the 
Anderston localization depends on the average that is 
used to calculate Pq(u>); using the geometrical average 
(q = 0), p q (io) vanishes for a certain value of A, while us- 
ing the arithmetical average (q = 1), the localization can- 
not be detected. Fig. [2]presents the phase diagram of the 
ground state for the Anderson- Falikov-Kimball model for 
U = 3.6 at t 2 = (solid line), t 2 — 0.4(dash-dot line) and 
t 2 = 1 (dash-dot-dot line) using the geometrical average 
(black lines) and the arithmetic average (brown lines). It 
is obtained observing the behavior of the average local 
density of states for different values of w as the disorder 
strength A is increased. We can identify three regions 
with respect to the states at the Fermi level: extended 
gapless phase (continuous spectrum), localized gapless 
phase (pure point spectrum) and gap phase @. The 
lines delimiting these phases are centered around w = 
for t 2 /t\ = 0, what does not happen for t 2 /t\ ^ 0. The 
asymmetry in the phase diagram becomes more evident, 
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FIG. 3: Averaged electron number n e versus chemical poten- 
tial at U = 3.6, for tt/t\ = and A = 0, 2.0, 4.0 and 6.0. 
(a)g = 1 and (b)q — 0. 
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FIG. 4: Averaged electron number n e versus chemical poten- 
tial at U = 3.6 using q = 0, for tj/tj / and A = 0, 2.0, 4.0 
and 6.0. (a) t* 2 /tl = 0.4 and (b) t%/t\ = 1.0. 



the higher the value of t^/t*. 

As already pointed out, if t% ^ 0, the chemical poten- 
tial ii = U/2 that was introduced does not affect the sys- 
tem in the half-filled band. When we look at the curves 
of Fig. [3] obtained for the relation between the chemical 
potential /i and the average number of electrons n, we see 
that all the curves obtained for i^A* =0 are centered at 
fi — U/2 and n e = 0.5, independent of the average that 
is considered. This means that the half-filled band of the 
system is always in a band center (u> = 0) as we vary the 
disorder A of the system. In Fig. 2]we present the curves 
corresponding to t 2 ^ {t 2 = 0.4 and i 2 = 1.0) using the 
geometrical average. We can observe in all curves that 
for a half-filled band system (n e — 0.5) the value of /i 
does not correspond to U/2 anymore, independent of the 
average that is used and the disorder that is imposed. 

Fig. shows the chemical potential versus disorder 
parameter A for the half-filled band (n e = 0.5). We 
observe that for t\/t\ — the half- filled band is in a 
band center {to = 0) independent of disorder A. As we 
include the influence of next-nearest-neighbors t\/t\ ^ 
the half-filled band is dislocated to bands with u> > 
that correspond to values of \i always greater than U/2. 

The ground-state phase diagram for electrons in a half- 
filled band as a function of t%/t\ for U — 3.6 is shown 
in Fig. [5] Performing the same process for other values 
of U, we can easily obtain all the simulation points of 



FIG. 5: Chemical potential versus disorder parameter A for 
the half-filled band (n e = 0.5). The results are obtained for 
t* 2 /tl = 0, 0.4 and 1.0. 



Fig. [7] In this figure we present a complete ground-state 
phase diagram for electrons in a half-filled band for three 
different values of for t^/t*, namely t\/t\ = 0, 0.4 and 
1.0. The presence of t 2 increases the extended phase. 
The dashed line represents the values of A and U where 
the energy gap goes through the half-filled band. Note 
that for ti different from zero and U large the extended 
phase does not disappear because the energy gap does 
not correspond to the half-filled band. 
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FIG. 6: Ground-state phase diagram for electrons in a half- 
filled band as a function of t%/t\ for U = 3.6 . Dots are 
determined from the numerical solution of the DMFT equa- 
tions. Lines are guide to the eye. 
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FIG. 7: Ground-state phase diagram for electrons in a half- 
filled band as a function of U for t^/t* = 0, 0.4 and 1.0. Dots 
are determined from the numerical solution of the DMFT 
equations. Lines are guide to the eye. 



FIG. 9: Ground-state phase diagram for electrons in a half- 
filled band versus a for U/t* = 0. Dots are determined from 
the numerical solution of the DMFT equations. Lines are 
guide to the eye. 



IV. LONG-RANGE CORRELATED DISORDER 

Recently results [l^] have shown that localization prop- 
erties are modified if correlations are introduced in the 
disorder distribution. For example, long-range correlated 
disorder favors the derealization, giving in ID systems 
a range of energies with extended eigenstates. The one- 
dimensional Anderson model with long-range correlated 
diagonal disorder displays a phase of extended electronic 
states where the on-site energy disorder distribution is 
described by a power-law spectral density. 

In what follows, we will study the effect due to 




FIG. 8: Ground-state phase diagram for electrons in a half- 
filled band by using long-range correlated disorder for a — 0, 
2.0 and 2.5. Dots are determined from the numerical solution 
of the DMFT equations. Lines are guide to the eye. 



long-range correlated disorder in the Anderson-Falicov- 
Kimball model on the Bcthc lattice using the DMFT ap- 
proach. The energy e; will not be considered as a random 
independent variable. Following an approach based on 
discrete Fourier transforms to construct long-range cor- 
related sequences, the on-site energies can be generated 
by the expression [ll| 
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where N is the number of sites and <f>k are N/2 indepen- 
dent random phases uniformly distributed in the interval 
[0,271"]. The uncorrelated disorder is recovered for a = 0. 

We will normalize the energy sequence to have the 
mean value < >= 0. The standard deviation is de- 
fined as a = a/< e| > - < e 4 > 2 = 6A/V3. For a = 0, 
we obtain a gaussian distribution for the . In this sense, 
the values of the Coulomb repulsion U and the disorder 
A that we obtain when po(w) = are not the same as 
obtained with the uniform distribution. 

In the case U = 0, we find in the literature that for the 
cubic network the transition from the extended states to 
the localized ones occurs for A c = 5.5 for the uniform 
distribution and A c = 7.0 for the gaussian distribution 
[22[ [23j . For the same value of U, we found in our sim- 
ulations A c = 5.7 and 6.8, respectively, for the uniform 
and the gaussian distributions. 

In Fig. [51 we present the ground-state phase diagram 
for electrons in a half-filled band by using long-range cor- 
related disorder and the geometric average. We see that 
for different values of a, the diagram that is obtained is 
similar to the one that corresponds to uncorrelated dis- 
order. 
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The difference between the diagrams can be quantified 
for different values of the exponent a when we look at 
the value of A when U — 0. The result is presented 
in Fig [3J We see that the value of A diminishes as the 
value of a is increased and tends to a constant value for 
a > 4. The decrease of A means that the correlation 
reduces the extended phase, different to what happens in 
low dimension [19(. 

V. CONCLUSIONS 

In the present paper, we studied the solutions of the 
Anderson-Falicov-Kimball model with nearest-neighbor 
(NNs) and next-nearest-neighbor (NNN) hopping. We 
found the averaged LDOS calculated using the arithmetic 
and geometric mean within dynamical mean field theory. 
We showed that the inclusion of the NNN with t* 2 jt\ ^ 
moves the half-filled band from the band center (uj = 0) 
to bands with lu > 0, independent of the disorder A 
that is considered. Besides that the chemical potential 
fi changes from U/2 for t^/t* = to greater values for 

We also analyzed the behavior of the system when 
long-range correlated disorder is considered using both 



the geometric and arithmetic means. We showed that 
the inclusion of this kind of disorder does not present 
results different from the ones already obtained for an 
uncorrelated disorder system At the other hand, we 
showed that decreasing the correlated disorder reduces 
the extended phase for U — 0, different to what happens 
in low dimension. 

A previous work about phase separation in the 
particle-hole assymetric Hubbard model involving NNN 
hoppings suggests that the inclusion of the temperature 
in determining the (U, T, fx) phase diagram is crucial for 
the study of the Anderson localization in a disordered 
electron system described by the Falicov-Kimball model. 
We will analyze these effects in future works. 
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